Entanglement between a qubit and the environment in the spin-boson model 
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The quantitative description of the quantum entanglement between a qubit and its environment is 
considered. Specifically, for the ground state of the spin-boson model, the entropy of entanglement 
of the spin is calculated as a function of a, the strength of the ohmic coupling to the environment, 
and e, the level asymmetry. This is done by a numerical renormalization group treatment of the 
related anisotropic Kondo model. For e = 0, the entanglement increases monotonically with a, until 
it becomes maximal for a — > For fixed e > 0, the entanglement is a maximum as a function of 
Q for a value, a — au < 1. 



Due to the promise of quantum computation there 
is currently considerable interest in the relationship be- 
tween entanglement, decoherence, entropy, and measure- 
ment. Motivated by quantum information theory several 
authors have recently investigated entanglement in quan- 
tum many-body systems 0,013,01 . It is often stated that 
decoherence or a measurement causes a system to become 
entangled with its environment. The purpose of this pa- 
per is to make these ideas quantitative by a study of 
the simplest possible model, the spin-boson model 0, 13 ■ 
This describes a qubit (two-level system) interacting with 
an infinite collection of harmonic oscillators that model 
the environment responsible for decoherence and dissi- 
pation. Specifically, we show how the entanglement be- 
tween a superposition state of the qubit and the envi- 
ronment changes as the coupling between the qubit and 
environment increases. One interesting result is that we 
find that the qubit becomes maximally entangled with 
the environment when the coupling a approaches a par- 
ticular finite value (a — > 1~). Furthermore, at this value 
the model undergoes a quantum phase transition, which 
is consistent with recent observations that often entan- 
glement is largest near quantum critical points 0,0,05 0- 

The spin-boson model. The Hamiltonian is 0,0] 

HsB = ^AcTj, 4- ^ea^ + ^ Wi(a|ai -t- i) 



(1) 



where A is the bare tunneling amplitude between the two 
quantum mechanical states t and i, £ is the level asym- 
metry (or bias), uji are the frequencies of the oscillators 
and Xi the strength with which they couple to the two 
quantum mechanical states. The effect of the oscillator 
bath is completely determined by the spectral function 



J{uj), defined belowQ. We will only consider the ohmic 
case, where it is has a linear dependence on frequency 

J{uj) = TT XfS{LU — LOi) — l-KClljJ 

% 

for <C Wc, and a is the dimensionless dissipation 
strength. The cutoff frequency, Wc ^ A. This model can 
describe the decoherence of Josephson junction qubits, 
such as those recently realized experimentally jTj, due to 
voltage fluctuations in the electronic circuit 0, and a can 
be expressed in terms of resistances and capacitances in 
the circuit and so this is an experimentally tunable pa- 
rameter. Recent results show it is possible to construct 
devices, with a <C I, the regime required for quantum 
computation. However, when modeling measurements 
one has a ^ I. 

The dynamical properties of the model have been ex- 
tensively studied. In particular, suppose the spin (qubit) 
is initially in a pure state which is a product state of 
up spin and the environment state, then the coherent 
Rabi oscillations that would be observed in the absence 
of coupling to the environment are modified as follows. 
One finds distinct behaviour for < a < 1/2 (damped 
coherent oscillations), 1/2 < a < 1 (exponential decay), 
and 1 < a (localization; i.e., the spin remains in the up 

state) HllEllS- 

Entropy of entanglement. We now consider a quanti- 
tative description of the entanglement of the qubit with 
the environment. A good entanglement measure for a 
pure state is the von Neumann entropy or entropy of 
entanglement 0, E| 



(2) 



where p is the reduced density matrix of the qubit. This 
is a two by two matrix given by 



\ a=x,y,z / 



(3) 
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where (tTa) denotes the expectation value in the state of 
interest. In this case Eq. ||2Jl reduces to 
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E{p) = -p+log2P+ -p- logaP-, 



(4) 
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where p± are the eigenvalues of the density matrix, 

p± = i(l±|(a)|). (5) 

For e = the only nonzero value of (aa) is {<Jx)- At 
T = it is given by 



dA ' 



(6) 



where Eq is the ground-state energy of the Hamiltonian 
(Q, and use has been made of the Feynman-Hellmann 
theorem. That the other values are zero can be seen by 
symmetry as follows. In general the Hamiltonian is in- 
variant under the reflection in spin space, ay ^<^y 
Hence, all eigenstates must have a definite parity under 
this transformation. Thus, {ay) — —{ay) for all states 
and so {ay) = at any temperature. For e — Q the 
Hamiltonian is also invariant under the joint transforma- 
tion — cTz and Qi — > ~ai and so {az) = at any 
temperature, provided there is no symmetry breaking. 

The challenge is now to evaluate the ground-state ex- 
pectation values {ax) and {az)- For a > 1/2 and par- 
ticularly for a ~ 1 this is a highly nontrivial problem 
because in this regime nonperturbative effects become 
importantly 13 . However, we show how these expectation 
values can be evaluated using the numerical renormaliza- 
tion group (NRG) applied to the equivalent anisotropic 
Kondo model. 

Anisotropic Kondo Model. The above model is equiva- 
lent to the anisotropic Kondo model (AKM), defined by 



k,a 

J\ 



^(cJ^coT - cl^coijS"" + gpBhSz- 



(7) 



The first term represents a free electron conduction band. 
We use a flat density of states po = 1/2-Do per spin, 
with 2Z?o the bandwidth. The second and third terms 
represent the transverse and longitudinal parts of the 
exchange interaction between a S" = 1/2 impurity and 
the local conduction-electron spin density, and the last 
term represents a Zeeman term for a magnetic field cou- 
pling only to the impurity spin. The correspondence be- 
tween H and Hsb, established via bosonization fl^ [l5l | . 
implies e = gpsh, ^ = poJj_ and a = (1 -I- 25/7r)^, 
where tan 6 = — 7rpo>/||/4. 5 is the phase shift for scat- 
tering of electrons from a potential J||/4 [sl. ITol. IT^ . We 
choose UJc = 2_Do so that A ~ J±. This equivalence 
has been used extensively to make predictions about the 
dynamics [lO| and thermodynamics [l5l | of the ohmic spin- 
boson model. The relevant low energy scale for the ther- 
modynamics is the Kondo scale Tk{ J± , J|| ) which is iden- 
tified with the renormalized tunneling amplitude, A^, of 
the spin-boson model. 



Ar 

UJc 



We restrict ourselves in this paper to the longitudinal 
sector of the AKM, i.e., to J± < \J\\\, where the sim- 
ple parameter correspondence between the models given 
above remains valid to lowest order in A/ujc = PqJ±- For 
larger values of A/ujc, a will acquire a renormalization 
due to finite A = Jj^, as indicated by the scaling analy- 
sis of the AKM in Rcfs. 13, This renormalization, 
however, is important mainly for the transverse sector of 
the AKM, J± > |J|||, which we do not consider in this 
paper. 

We turn now to the evaluation of {a^)- The equivalence 
between models ensures that the AKM has (to within an 
additive constant) the same ground-state energy Eq as 
that of the spin-boson model. At T = 0, we therefore 
find, in analogy to Eq. JHJl applied to the AKM with 
A = J_L, that 



{<^x) = { 4 iCcx^" + H.c.) 



(9) 



(8) 



i.e., {ax) can be obtained from a local static correla- 
tion function. Another way of seeing that this relation 
is valid, is to note that the unitary transformation in 
bosonization which transforms H into Hsb also trans- 
forms (cj ^co^iS~ + H.c.) into ax of the spin-boson model 
(details of this mapping can be found in Ref. 0| and in 
greater detail in Appendix A of Ref. ITsil . The same uni- 
tary transformation on the AKM transforms (5'^) into 
tTz/2 of the spin-boson model. The latter can therefore 
be calculated directly within the AKM as a thermody- 
namic average {2Sz). 

Method. The above local correlation function can be 
calculated from Wilson's NRG method 16] which has 
been shown to give very reliable results for quantum 
impurity models such as the AKM|0|. The approach 
used here allows in addition the calculation of local dy- 
namical quantities, such as the dynamical susceptibility 
(( cl^co^iS--cl^CQ^^S+)) m. In outline (see Ref. [H 
for the details), the procedure consists of introducing a 
logarithmic mesh of k points A:„ = A~",A > 1 for the 
conduction band and performing a unitary transforma- 
tion of the Cfccr such that /oo- = J2k '^ka is the first op- 
erator in a new basis, fncnn = 0,1,..., which tridiag- 
onalizes He — X^fe^t ^kfj.cl^Ck^ in k space. The Hamil- 
tonian |(7J) with the new discretizcd form of the kinetic 
energy is now diagonalized by the following iterative pro- 
cess: (a) one defines a sequence of finite size Hamil- 

E^EL"o enA-"/2(/t^^^/„^ + H.c.) + 

/oVot5+) + 4^(^t/oT - flfoi)S'' for 
iV > and ^„ 1 for n 1 jl^; (b) the sequence of 
Hamiltonians Hn for TV = 0, 1, ... is iteratively diagonal- 
ized within a product basis of, typically, up to 1200 states 
for each iteration, up to a maximum value N = Nm. 
This gives the excitations and many body eigenstates 
at a corresponding set of energy scales ujn defined by 

N — 1 

the lowest scale ujn — A ~ in Hpf. The matrix el- 
ements {m\Ox.z\n)N for the operators Ox = ^ coals'" 
and Oz — Sz, required to calculate (ax) and (cr^), are 



N 



tonians H 
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FIG. 1: The dependence of (i) the ground-state expectation 
value {(Jx) as a function of the dimensionless couphng a to the 
environment for e = 0, and (ii) the entanglement entropy E of 
a qubit ohmically coupled to an environment as a function of a 
for £ = 0. The different curves correspond to different values 
of the ratio of the bare tunneling amplitude A to the cutoff 
frequency of the boson bath ujc- Note that as a — > 1~ the 
qubit becomes maximally entangled with the environment. 



also calculated iteratively. The choice of Nm depends on 
the Kondo scale Tk = and hence on a, but for given a 
(i.e. for given J±, J\\) should be large enough such that 
wtVto ^ A discretization parameter A = 1.5 was 

used throughout and we checked that the above expec- 
tation values remained unchanged on further increasing 
Nm- This suggested that our approximation of using a 
finite Nm to calculate the thermodynamic expectation 
values of the infinite system is a very good one. 

Symmetric case. For e = and a < 1 only (ax) is 
nonvanishing. We show this in Fig. Q at T = ver- 
sus dimensionless dissipation strength a in the range 
< a < 1, and for several values of the dimension- 
less tunneling amplitude A/ojc- The limiting noninter- 
acting value (ax) ^ 1 is recovered as a ^ 0. In the 
limit A /ujc ^ it vanishes at the quantum critical point 
of the spin-boson model a — 1 where A^ — > 0. For 
any finite fixed A/lUc, however, our use of the AKM im- 
plies that the critical behaviour occurs at ac > 1 with 
Qfc ^ 1 as J_L —> (specifically, this critical behaviour 
occurs at the ferromagnetic-antiferromagnetic boundary 
J± ~ Figure |^ also shows the T = entropy 

of entanglement of the qubit. The entropy vanishes as 
a — > and a ppr oaches its maximum value as a — > 1^ 
(see also Ref. [l^ for weak dissipation results for E). For 
q; > 1, we are in the ferromagnetic sector of the AKM 
where (cr^) — 1, {cFx) = 0, and the reduced density matrix 
eigenvalues p± = 0, 1 giving E = 0, i.e., E{a) drops dis- 
continuously at the quantum critical point a = 1 |2Q|] . It 




a 



FIG. 2: The dependence of the entanglement entropy of the 
ground state on the coupling to the environment a and the 
level asymmetry e for A/uJc — 0.04. Note that for e > 0, the 
entanglement is a maximum at a = um < 1- 



is interesting that for a spin qubit coupled to two bosonic 
baths it is possible to remain in the delocalized phase 
(i.e., (ctz) = 0) for all dissipation strengths j2j|. Finally, 
we note that the entropy of entanglement is quite differ- 
ent from the thermodynamic entropy of the boundary (or 
impurity spin entropy). The latter is usually defined as 
S{a) — S{a = 0) where S{a) is the total thermodynamic 
entropy of the system . The impurity spin entropy is 
zero for a < 1 because the ground state of the AKM is a 
spin singlet for J|| > 0. 

Asymmetric case. For e > 0, (ctz) acquires a finite 
value analogous to the magnetization (Sz) in a local mag- 
netic field gfish = £ in the AKM. The entanglement en- 
tropy E now depends on |(ct)| = (Kctk)]^ -I- [(cz)!)^^^ via 
Eq. © and is shown in Fig. Q . The behaviour of E as 
a function of a and e is understood from the behaviour 
of (ax) and (ctz) shown in Fig. (jSJ. In particular, we now 
find that for arbitrary small e, the entanglement entropy 
first increases with increasing a before reaching a maxi- 
mum value at a = um < 1 and then decreasing as a — >■ 1. 
This behaviour arises from the competition between the 
a dependence of (ax) and (cr^) in Fig. Q. Whereas (ax) 
continues to decrease monotonically with increasing a (as 
for e = 0), it is seen that {(Tz) increases monotonically 
with increasing a with (az) ^ 1 as a — > 1. The condi- 
tion for full polarization, (cr^) « 1, is e 3> A^. For any 
e > 0, this condition is always satisfied since A^ — > as 
a — > 1. It follows that |((t)| has a minimum as a function 
of a and that the entanglement entropy, for e > 0, has 
a maximum at a — am < 1 before decreasing again as 
a ^ 1. 

Finally, we suggest several directions for future work, 
(i) This work focused solely on static properties of the 
spin-boson model. It would be interesting to consider 
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FIG. 3: The dependence of the ground state expectation 
values (ax), (ctz) on a and e for A/iOc = 0.04. For a — 
the noninteracting values (a^) = A/y/e^ + and {(Tz) — 
e/\/e^ + A^ axe recovered for all values of e/A. 

dynamics, for example, the longitudinal and the trans- 
verse dynamical susceptibilities, and hence extract the 
decoherence and relaxation rates for an ohmically cou- 
pled qubit. In addition, it is interesting to ask how the 
entanglement varies vifith time if the initial state has no 
entanglement of the qubit and the environment. 

(ii) The AKM is integrable by the Bethe ansatzp^j. 



The AKM can also be related to a free boson field theory 
with a boundary sine-Gordon term|^ which is also 
integrable by the Bethe ansatz. Exact expressions can 
be obtained for the free energy. It involves solving a set 
of thermodynamic Bethe ansatz (TBA) equations. At 
T = the impurity ground state energy is going to be 
related to Tk- The real problem is getting results for 
arbitrary anisotropics (dissipation strengths) 24] . 

(iii) Recently it was shown^ that if in a quantum crit- 
ical system one calculates the entropy of entanglement of 
a subsystem of size L with the rest of the system this 
equals the geometric entropy previously calculated for 
the corresponding conformal field theory (motivated by 
questions concerning black hole thermodynamics!) [2^. It 
would be interesting to perform similar calculations for 
the relevant boundary field theory. 

(iv) The NRG can also be used to reliably cal- 
culate properties of the spin-boson model at nonzero 
temperature :2Qj. However, calculating the entanglement 
at nonzero temperature is an open problem because it 
involves a mixed state and it is not practical to evaluate 
the measure of entanglement that has been proposed for 
such states 
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